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The superfluidity of Bose-Einstein condensates in optical lattices are investigated. Apart from the 
usual Landau instability which occurs when a BEC flows faster than the speed of sound, the BEC 
can also suffer a dynamical instability, resulting in period doubling and other sorts of symmetry 
breaking of the system. Such an instability may play a crucial role in the dissipative motion of a 
trapped BEC in an optical lattice recently observed fll. 



PACS numbers: 67.40.Db, 03.75.B, 05.30.Jp, 32.80.Pj 

Optical lattices have long been used to manipulate 
ultra-cold atoms, with applications ranging from beam 
splitters ||^, accelerators [^ ^ to lithography Q. There 
are now growing interests in replacing the cold atoms 
with Bose-Einstein condensates (BECs) of alkali atoms 
[g|, 0, H to explore the effects of coherence, atomic inter- 
action, and superfluidity, with important applications in 
atom lasers |9|] and high-precision interferometry |10|| . 

In this Letter, we investigate the superfluidity and in- 
stabilities of BECs in optical lattices. In free space, the 
superflow of a uniform BEC is represented by a plane 
wave, which has Landau instability when it travels faster 
than the sound. In an optical lattice, the natural ob- 
jects of concern are the BEC Bloch waves, whose am- 
plitudes are modulated with the same periodicity as the 
lattice. In the central region of the Brillouin zone of the 
lowest band, they are found to be local energy minima 
against all sorts of perturbations and thus represent the 
superflows of the BEC in optical lattices. For sufRciently 
strong repulsive interactions between the atoms, the su- 
perfluidity region can extend over the entire Brillouin 
zone; but for weaker interaction. Landau instability can 
occur in the outer regions of the zone. Many of the Bloch 
states with Landau instability can even be dynamically 
unstable in that small initial disturbances around them 
grow exponentially in time, resulting in period doubling 
and other forms of spontaneous breaking of the periodic- 
ity of the system. This dynamical instability is unique to 
BEC Bloch waves and not present in BEC plane waves 
in free space. We map out the dangerous zones of the 
dynamical instability, characterize the growth rates, and 
discuss the experimental consequences. 

We consider the situation of a one dimensional optical 
lattice in which the motion in the perpendicular direc- 
tions are confined |Q] or can be disregarded B . We treat 
the atomic interaction with the mean field theory, and 
obtain the grand canonical Hamiltonian 
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where all the variables are scaled to be dimensionless by 
the system's basic parameters, the atomic mass m, the 
wave number fc^ of the two laser lights that generate the 
optical lattice, and the average density no of the BEC. 
The chemical potential fi and the strength of the periodic 
potential v are in units of 
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the wave function in 
and t in units of ™^ . 

where Oc > is the 



units of ^yno, x in units of 
The coupling constant c = 
s-wave scattering length. 

Hamiltonian (|l|) is extremized by states in the form of 
Bloch waves, (j)k{x) — e^'^^ipkix), where ipk{x) is of the 
period of the optical lattice and can be expanded as a 
Fourier series. To find the numerical solution of a Bloch 
state in the lowest band, ipkix), we truncate the series 
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up to A^th term (pk{x) — X^-jv*^™^""^ (^^ ^^^"^ ^ ^ 
10). The numerical solution is obtained by varying {um} 
so that the wave function ^pk{x) minimizes the system's 
total energy. The accuracy is checked by substituting the 
solutions into the Gross-Pitaevskii equation, ~2^J2'^"'" 
V cos X (/)-|~c|(/)p0 = /i0 , which is obtained by the variation 
of Hamiltonian (^ . 

To determine the superfluidity of these Bloch states, 
we need to find out if they remain energy minima against 
perturbations which break the periodicity. These pertur- 
bations can be decomposed into different modes labeled 
by 9, 



5Lpk{x, q) = Uk{x, q)e"''' + vl{x, q)e '«^ 



(2) 



where q ranges between —1/2 and 1/2 and the pertur- 
bation functions Uk and Vk are of periodicity of 27r in x. 
Since the system is periodic, the quadratic form of the en- 
ergy deviation from the Bloch state 4>k is block diagonal 
in 5, with each block given by 



where 
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with 



v=0.01 



v=0.05 



v=0.10 
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If Mk{q) is positive definite for all —1/2 < q < 1/2, the 
Bloch wave (jik is a local minimuni. Otherwise, SEk can 
be negative for some q, and the Bloch wave is a saddle 
point. 

We first consider the special case t; = 0, BEG in free 
space, where the Bloch state (jjk becomes a plane wave 



^ikx 



The operator Mk{q) becomes a 2 x 2 matrix 



Mk{q)^ 



q^/2 + kq - 



q" /2 — kq + c 



whose eigenvalues are found easily as 



A± = 



;± y/k^q 
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Since A+ is always positive, M]^{q) fails to be positive def- 
inite only when A_ < 0, or equivalently, \k\> ^Jq^ /A + c. 
It immediately follows that the BEG flow e^^^ becomes 
a saddle point when the flow speed exceeds the sound 
speed, |fc| > ^/c. This is exactly the Landau condition 
for the breakdown of superfluidity pl| , which has re- 
cently been confirmed experimentally on BEG llj] . 

The stability phase diagrams for BEG Bloch waves are 
shown in the panels of Fig.l, where different values of v 
and c are considered. The results have reflection sym- 
metry in k and g, so we only show the parameter region, 
< A: < 1/2 and < q < 1/2. In the shaded area (light 
or dark) of each panel, the matrix Mk{q) has negative 
eigenvalues, and the corresponding Bloch states <f>k are 
saddle points and have Landau instability. For those val- 
ues of k outside the shaded area, the Bloch states are local 
energy minima and represent super-flows. The super- flow 
region expands with increasing atomic interaction c, and 
occupies the entire Brillouin zone for sufficiently large c. 
On the other hand, the lattice potential strength v does 
not affect the super-flow region very much as we see in 
each row. The phase boundaries for v <^1 are well repro- 
duced from the analytical expression k = y'q^jA + c for 
V = 0, which is plotted as triangles in the first column. 

A saddle-point Bloch state (pk can still be dynamically 
stable in that small deviations from it remain small in 
the course of time evolution if no external persistent per- 
turbations are present. This is the case for all Bloch 
states either in the absence of atomic interactions or pe- 
riodic potentials. When both factors are present, many 
of the saddle-point Bloch-states become dynamically un- 
stable against certain perturbation modes q, shown as the 
dark-shaded regions in Fig.|l|. These results are obtained 
from the linear stability analysis of the Gross-Pitaevskii 
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FIG. 1: Stability phase diagram of BEG Bloch states in op- 
tical lattices, k is the wave number of BEG Bloch states; 
q denotes the wave number of perturbation modes. In the 
shaded (light or dark) area, the perturbation mode has neg- 
ative excitation energy; in the dark shaded area, the mode 
grows or decays exponentially in time. The triangles in (a.l- 
a.4) represent the boundary, g^/4 -I- c = k^, of saddle point 
regions at u = 0. The solid dots in the first column are from 
the analytical results of Eq.(h2|). The circles in (b.l) and (c.l) 
are based on the analytical expression (nsl) . The dashed lines 
indicate the most unstable modes for each Bloch state k. 



which governs the dynamics of the system. A Bloch state 
(pk is a stationary solution of this equation, depending on 
time only through the phase factor e~''^'. Writing the 
deviation in the similar form of Eq.(S), and expanding 
the above equation to first order, we find 



d 



dt V Wfc 
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The dynamical stability of the Bloch state (pk is deter- 
mined by the eigenvalues ek{q) of the matrix aMk{q). If 
they are real for all —1/2 < q < 1/2, the state is dynam- 
ically stable; otherwise it is dynamically unstable. 

Before discussing our detailed results on the dynamical 
instability, we pause here to make some general remarks. 



(i) When all the eigenvalues ek{q) are real, the motions 
around the Bloch state 0^ are oscillations, which can be 
quantized to yield the phonon excitations |0, O. The 
traditional Bogoliubov approach yields the same matrix, 
aMk{q), for the phonon spectrum. However, the bosonic 
commutation relation for the phonon operators imposes 
the skewed normalization condition X'^crX = 1, which 
selects only half of all the modes. The other half, satis- 
fying X''(tX = — 1, will be called antiphonon modes for 
ease of reference, but they really do not represent phys- 
ical degrees of freedom independent of the phonons. (ii) 
When the Bloch state is a local minimum {Mk{q) posi- 
tive definite), the dynamical eigenvalues ek{q) oiaMkiq) 
are all real and the phonon branch of the spectrum is 
positive. The key to the proof is to notice that 



ek{q)X^aX = Xnh{q)X 



(10) 



for an eigenvector X of aMk (q) ■ Because the right hand 
side is positive and X^aX is real, efc(q) is real and has 
the same sign as X^aX. The physical meaning of this 
theorem is that, when it is a local minimum, the Bloch 
state (f>k is dynamically stable and its phonon excitations 
are not energetically favored, (iii) Because the matrix 
aMk{q) is non-hermitian and real (when expressed in 
the momentum representation) , complex eigenvalues can 
only appear in conjugate pairs, corresponding to modes 
growing or decaying exponentially at rates given by the 
imaginary part of the eigenvalues. Because both the 
quadratic forms in Eq.(|l^ are real, they must vanish 
when Skiq) is complex. It is then impossible to enforce 
the normalization condition X'^aX ~ 1, corresponding 
to the fact that such modes cannot be quantized. 

We now present our detailed results on the dynamical 
stability. Again, we first look at the case ti = 0, where 
the eigenvalues of a Alt (q) are 



£±(g) = kq± Vrc + 9V4- 



(11) 



These eigenvalues are always real; the BEC flows in free 
space are always dynamically stable. When v ^ 0, the 
situation is totally different: the eigenvalues ek{q) of 
aMk {q) can be complex and Bloch states can be dynam- 
ically unstable. The dark-shaded areas in Fig. ^ are the 
places where these Sk{q) are complex. 

In the first column of Fig.|^, where u <C 1, the dark- 
shaded areas are like broadened curves. These curves are 
the solutions of e+{q — 1) = £-(?), 



k = ^q^c + qyA + ^{q - l)^c + (q - 1)4/4, (12) 

which are plotted as solid dots in Fig.^. This is the 
resonant condition for a phonon mode to couple with an 
antiphonon mode by first order Bragg scattering. The 
resonance is necessary because the complex eigenvalues 
can appear only in pairs, and they must come from a pair 
of real degenerated eigenvalues. Resonances within the 



phonon spectrum or within the antiphonon spectrum do 
not give rise to dynamical instability; they only generate 
gaps in the spectra. Somehow, in order to produce a 
mode with zero normalization, one must couple a pair of 
modes with opposite normalizations. 

In the first row of Fig. |[ we have another extreme case 
c <C w. The open circles along the left edges of these two 
dark-shaded areas are given by 

Eiik + q)- Eiik) = Ei{k) ^Ei{k- q) (13) 

where Ei{k) is the lowest Bloch band of Hq — —5^^ + 
wcosx. When c = 0, this periodic system is linear; the 
excitation spectrum just corresponds to transitions from 
the the condensate of energy Ei (fc) to other Bloch states 
of energy En{k -l- 5), or vise versa. The above equation 
is just the resonant condition between such excitations 
in the lowest band {n = 1). Alternatively, this condi- 
tion may be viewed as the energy and momentum con- 
servation for two particles in the condensate to interact 
and decay into two different Bloch states Ei {k + q) and 
Ei{k-q). 

One common feature of all the diagrams in Fig.|l| is 
that there is a critical Bloch wave number kd beyond 
which the Bloch states (t)k are dynamically unstable. The 
onset instability at kd always corresponds to q = 1/2. In 
other words, if we drive the Bloch state 0^ from A; = to 
k = 1/2 the first unstable mode appearing is always q — 
±1/2, which represents period doubling. Only for kkd 
can longer wavelength instabilities occur. The growth 
of these unstable modes drives the system far away from 
the Bloch state and spontaneously breaks the translation 
symmetry of the system. The critical value of the Bloch 
wave number for the case oiv <^ 1 is found to be kd = {c-\- 



1/16)^ 



into 



Eq.(|T|). 



In the other 



by substituting q — 2 
extreme case, c ^ u, the same substitution in Eq.(p^) 
yields kd — 1/4 with the help of periodicity of the band 
energy. Based on these results and the diagrams in Fig.[l|, 
we find that kd> \- 

The dynamical instability discovered in this work 
should be observable in experiments. We have mapped 
out the dangerous zones of dynamical instability, which 
give us a good sense of where to look for unstable Bloch 
states and modes of instability. In Fig.||, the rate of 
growth r for the most prominent mode (dashed lines in 
Fig. P of each Bloch state k is plotted in Fig.g. The 



physical unit of the growth rate is 



Ahk, 



which is 4.0 



per microsecond for sodium and 0.16 per microsecond 
for rubidium. Since the lifetime of BECs can be up to 
the order of seconds |15 , these growth rates in Fig. Hare 
significant. It is possible to directly observe the change 
of periodicity of the BEC due to the dynamical instabil- 
ity by monitoring the Bragg scattering of a probing laser 
light by the BEC cloud ||l3|. This dynamical instability 
can also cause the disruption of Bloch oscillations ^\ . 
In a recent experiment, the superfluidity and instabil- 
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FIG. 2: Growth rates of the most unstable modes of each 
Bloch state <j)k- The erratic behavior of some curves around 
k — 1/2 is due to the difficulty in finding the accurate Bloch 
waves (/>fe in this region when c > v |17[. 



ities of a BEC in an optical lattice was studied using a 
cigar-shaped (one dimensional) magnetic trap R]. Af- 
ter the BEC was prepared in the trap in the presence 
of the optical lattice, the trap was suddenly shifted by 
Aa; along the longitudinal direction. This is equivalent 
to displacing the whole BEC off the center of the har- 
monic trap then releasing it. The subsequent oscillations 
of the BEC are non-damped if the initial displacement 
is small, but become dissipative if Ax is over a criti- 
cal value Axc- This qualitatively agrees with our stabil- 
ity diagrams, because larger Ax implies larger velocity 
and therefore larger k. The dissipative behavior was ex- 
plained as a manifestation of the Landau instability, but 
the dynamical instability discussed here is likely to play a 
crucial role in our view. First, the experiment has v '^ 0.2 
and c ^ 0.02, where the dynamical instability is rampant 
according to Fig. 2. Second, Axc increases with decreas- 
ing lattice potential v, which is in accordance with the 
trend of the growth rate as a function of v and k shown in 
the figure. Third, there is no dissipation when the BEC 



density is low, where Landau instability should be very 
strong but dynamical instability be very weak according 
to our Fig.l. However, more detailed analyses are needed 
to take into account of the effects of inhomogeneity and 
thermal cloud before a quantitative comparison with the 
experiment. 
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